Ferroelectricity of Ice Nanotubes inside Carbon Nanotubes 
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We report that ice nanotubes with odd number of side faces inside carbon nanotubes exhibit 
spontaneous electric polarization along its axis direction by using molecular dynamics simulations. 
The mechanism of this nanoscale quasi-one-dimensional ferroelectricity is due to low dimensional 
confinement and the orientational order of hydrogen bonds. These ferroelectric fiber structural 
materials are different from traditional perovskite structural bulk materials. 
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PACS numbers: 77.84.-s, 85.35.Kt, 61.46.Fg, 68.08.De 

Whether ferroelectric ice exists is a question that has 
long fascinated researchers [J H, H, 0]. Although one wa- 
ter molecule has a significant dipole moment, normal 
hexagonal ice Ih does not possess ferroelectricity and it 
is very difficult to get an ice phase even with a weak 
whole polarization. Several year ago, the ferroelectricity 
of ice films grown on ultra-clean platinum was detected, 
although only a small proportion (0.2%) of the molecules 
are aligned [a, @. Most recently, Fukazawa et al. have 
succeeded in making Ice XI in the laboratory and sug- 
gests existence of ferroelectric ice in the universeQ. 

Confinement of matter on the nanometer scale can in- 
duce phase transitions not seen in bulk systems 0, Q. 
In biology as well as in natural and synthetic materi- 
als, water is often tucked away in tiny crevices inside 
proteins or in porous materials. Therefore, water con- 
fined in nanoscale quasi-one-dimensional (Q1D) channels 
is of great interest to biology, geology, and materials sci- 
ence. Recent investigations have increased our under- 
standing of confined water, showing that, in nanoscopic 
proportions, many water properties differ drastically 
from those of bulk water [ToL ITU IT^ . In particular, 
an excellent model is water confined in carbon nan- 
otubes, which have been adopted in many previous stud- 
ies [IlQITllMITSEITilM 

Among these interesting investigations, Koga et al. in- 
vestigated water inside single-walled carbon nanotubes 
(SWCNTs) of diameter 1.1—1.4 nm and found that wa- 
ter forms ice nanotubes composed by a rolled square ice 
sheet 0, [H EH. Recent neutron scattering studies, 
combined with molecular dynamics (MD) simulations, re- 
vealed that water in SWCNTs of diameter 1.4 nm forms 
a core-shell structure, in which the square ice sheet de- 
scribed by Koga et al. is coupled with a chainlike config- 
uration at the center of the shell [27], [28j] . 

Recently, investigations of nanoscale ferroelectric ma- 
terials are very active due to its potential application 
[29L I30L I3TI ] . Because of its particular structures, the ice 
nanotubes inside SWCNTs may exhibit ferroelectricity. 



In this paper, we report that ice nanotubes with odd 
number of side faces inside SWCNTs exhibit sponta- 
neous electric polarization along its axis direction by us- 
ing MD simulations. The origin of ferroelectricity of these 
ice nanotubes is due to the strong direction preferred 
hydrogen-bonds between water molecules and the Q1D 
confinement inside SWCNTs, which is different from tra- 
ditional ferroelectric bulk materials, such as displacive 
perovskite-based oxides. 

The water-water interaction is modelled by the TIP4P 
model j32j, which has been used in many simulations 
and given theoretical results in well consistent with ex- 
periments [13, IH HE HI. The Ewald summation 
is used to deal with the long range Coulomb interac- 
tion and the Lennard- Jones interaction between water 
molecules is cut off at 9.0 A. The interactions between 
water molecules and SWCNTs are contributed only by 
the Lennard-Jones potential between the oxygen and car- 
bon atoms 18], which can be written as 



N„ N c 



(1) 



i=l j=l 



where N a is the total number of oxygen atoms, N c is the 
total number of carbon atoms, and ry is the distance 
between the oxygen and carbon atoms. The parameters 
are taken as e oc = 0.1029 kcal mol -1 and a oc = 3.280 
A, obtained from Lorentz-Bcrthlot combining rules of 
foe = \Je QO x e cc and a oc = (a oa + cr cc )/2. The pa- 
rameters, e 00 — 0.1549 kcal mol -1 and a oa — 3.154 A 
are taken from the TIP4P model, e cc — 0.06835 kcal 
mol" 1 and <j cc = 3.407 A are taken from multiple-shell 
fullerences and multiwalled carbon nanotubes [3 [35J. To 
speed up MD simulations, we assumed that the carbon 
atoms on SWCNTs are distributed continuously, which 
is well known as the continuum model and has been used 
in a lot of systems 36] . Thus the potential between one 
oxygen atom and a SWCNT is written as 



U(r) — rid u{x)dY,, 
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(2) 



where r and x represent the distances of the oxygen atom 
to the axis of SWCNT and the surface element cE, re- 
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FIG. 1: (color online) Relaxed structures of ice nanotubes in- 
side SWCNTs. The red and gray balls denote oxygen and hy- 
drogen atoms, respectively. The ice nanotubes inside (15,0) 
and (17,0) SWCNTs have spontaneous electric polarizations 
along its axis direction. The polarization vectors are illus- 
trated by P z and arrows. 



spectively. rid is the mean surface density of atoms on 
SWCNT HI HI. The c " c bond lcn S th of SWCNTs is 
assumed to be 1.42 A. We checked the difference between 
the formula (1) and (2) is within 2% in our simulation 
conditions. 

The MD simulations are performed by a modified TIN- 
KER package [39( , in which the velocity- Verlet algorithm 
for integration of equation of motion and the Berendsen 
thermostat algorithm are used. The time step is set to 
be 0.5 fs. The motion of translation and rotation of the 
whole water molecules are eliminated during the MD sim- 
ulations to calculate temperature. The simulation box is 
taken as a 50Ax50AxL cubic box, where L is the length 
along the axis of SWCNT. The periodic boundary con- 
ditions along the axis are used. Then, the real axial- 
pressure (P zz ) is the original calculated value times 2500 

A /S, where S is the area of the cross section of SWCNT. 

Following Koga et al. [H, [3, [HI], we first simu- 
lated water inside SWCNTs at constant temperature 
(T) and P zz by using NPT assemble. We chose a high 
axial-pressure (P Z2 =200MPa) condition and four types 
of zigzag SWCNTs, which are (14,0), (15,0), (16,0), and 
(17,0) with diameters of 10.96 A, 11.74 A, 12.53 A, and 
13.31 A, respectively. During our MD simulations, the 
temperature was lowered stepwise staring from 450K to 
100K and then heated stepwise back. The MD simula- 
tion time at each temperature is 2ns (10ns or more near 
the phase transition point). Some relaxed structures of 
the ice nanotubes and the polarized directions are both 
shown in Fig.l. These results consistent with those of 
Koga et a/. [l3, [H, [IE Ell- As illustrated in Fig.l, the ice 
nanotubes have 4, 5, and 7 side faces corresponding to 
inside (14,0), (15,0), and (17,0) SWCNTs, respectively. 
The ice nanotube inside (16,0) has 6 side faces but not 
shown in Fig.l. From the viewpoint of application, we 
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FIG. 2: (color online)The potential energy per water molecule 
(Epot, in which the water-SWCNT interaction energy is ex- 
cluded) and polarization per layer (P z ) versus temperature. 
The left and right panels are of the ice nanotubes inside (15,0) 
and (17,0) SWCNTs, respectively. The red filled circles and 
blue unfilled circles indicate the heating and cooling processes, 
respectively. The statistic standard errors are shown by the 
error bars, and lines are to guide the eye. 



think that water or ice will probably be enclosed inside 
SWCNTs to produce Q1D devices. So, our later MD 
simulations were performed at constant T and volume 
(V or L), known as NVT assemble. The initial densities 
are taken by considering the results at high P zz . In our 
simulations, the L is taken as 109.2 A and the total num- 
ber of water molecules in the simulation box is 40 x n, 
where n is the number of side faces and refers 4, 5, 6 and 
7 corresponding to the four types of ice nanotubes inside 
(14,0), (15,0), (16,0), and (17,0) SWCNTs. The relaxed 
structures of the ice nanotubes and the polarized direc- 
tion are almost the same as those at high P zz shown in 
Fig.l. It is found that the ice nanotubes with odd num- 
ber of side faces are ferroelectric, while those with even 
number of side faces are antiferroelectric. The polariza- 
tions of the relaxed ice nanotubes with n=5 and 7 are 
0.351 eA and 0.347 eA per layer of water molecules. If 
we take a bundle of these ferroelectric "fibers" arranged 
in a dense packed hexagonal 2D lattice with a interspace 
distance of 3.4 A (the same distance as that in graphite), 
the polarized rates are 0.15 C m -2 and 0.12 C m~ 2 , re- 
spectively (the rate of BaTi03 is 0.26 C m~ 2 ). 

The potential energy (E pot ) and spontaneous electric 
polarization along axis (P z ) versus T of two ferroelectric 
ice nanotubes are plotted in Fig. 2. It is found that the 
ferroelectric-paraelectric of the ice nanotubes strongly 
depends on the solid-liquid phase transition. The fer- 
roelectric phase transition behavior of the ice nanotubes 
is much different from traditional ferroelectric bulk ma- 
terials such as BaTiOa, whose transition occurs at solid 
state phase and obeys the Curie- Weiss law. As an order 
parameter, the P z of an ice nanotube drops at the solid- 
liquid phase transition point consisting with its E pot -T 
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FIG. 3: (color online) Sketch of an unrolled 2D ice sheet, in 
which, the basis vectors (ai and 3%), the chiral vector (Ch), 
and the axis vector (T) are shown. Ch = nai + ma2 corre- 
sponds to a (n, m) ice nanotube. The "Bernal-Cowler-Pauling 
bulk ice rule" is that every water molecule serves as a double 
donor and a double acceptor of hydrogen bonds, and every 
water molecule is hydrogen-boned to exactly four nearest- 
neighbor molecules0, In the 2D ice sheet and ice nan- 
otube, the "ice rule" is still satisfied and the water network 
looks like a chessboard consisting of many square lattices. The 
square hydrogen-bond network is denoted in dashed lines and 
the details of the hydrogen bonds are shown at the center 
by four water molecules. The red arrow at a site denotes 
the dipole moment of a water molecule and is to direct its 
orientation. The example ice nanotube in this figure is (5, 1) . 



curve. This ferroelectric-paraelectric phase transition is 
a first-order transition coupled with violent structural 
change. The ferroelectricity is due to the special struc- 
tures of ice nanotubes and the orientations of water 
molecules, which details are discussed as follows. 

Just like for SWCNTs, we can analyze the structure of 
an ice nanotube by rolling a 2D ice sheet. For simplicity, 
we use a right-angle water molecular model to describe 
the 2D ice sheet. As illustrated in Fig. 3, a 2D ice sheet 
consists of many square lattices and an ice nanotube can 
be described by a pair of numbers (n, m). By using this 
description, the three ice nanotubes shown in Fig.l can 
be denoted as (4, 0), (5, 0), and (7, 0), respectively. A pair 
of (n, m) can describe the main frame of water network 
of an ice nanotube but fail to describe the orientations of 
water molecules. If we use the dipole moment of a water 
molecule p (0.45 eA in the TIP4P model and 0.38 eA in 
experiment) to direct its orientation, we can get another 
equivalent description of the "ice rule" in a 2D ice sheet 
and so in an ice nanotube: the projections of all diploe 
moments aligned along the basis vectors (d\ and 3,2 in 
Fig. 3) must be in the same direction. That is to say, 
there are only two types of orientation for two nearest 
water molecules, parallel or perpendicularly. 

To understand why the water molecules prefer the ori- 
entations shown in Fig.l, a very simple model is used 
to compare the energies of different isomers. We take 
the approximation: En ~ Elj + Ehy + Ep, here Elj is 
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FIG. 4: (color online)The upper panel: four types of orienta- 
tions of water molecules consisting one side face in a (n, 0) ice 
nanotube; The lower panel: unrolled 2D hydrogen-bond net- 
work of two example ice nanotubes, (4, 0) and (5, 0). The red 
arrows indicate the orientations of water molecules as show 
in Fig. 3. The potential energies of (a), (b), (c) and (d) are 
in increasing order, so the type (a) is preferred. A (4, 0) ice 
nanotube has 4 water molecules per layer with 2 upwards and 
2 downwards, so it has no polarization at every layer. Every 
layer of a (5, 0) ice nanotube has 5 water molecules with 3 
upwards and 2 downwards and has a polarization along the 
axis due to the one redundant upward water molecule. All 
the layers have the same directions of polarization along the 
axis since the "ice rule" forbids any opposite one. The "ice 
rule" accumulates the polarizations at every layer and gener- 
ates a whole spontaneous electric polarization P z along the 
axis. The blue boxes are to denote the unit cell. 



the Lennard- Jones potential, Ehy is the potential owned 
by hydrogen bonds, and the Ep is the potential due to 
dipole-dipole interaction which can be written as 



E f 



1 



47re ( 



■Pi - Khj -PiWij -Pj) 



(3) 



where eo is the permittivity of space, pi and pj are the 
dipole moments of water i and j, respectively, is the 
displacement from water i to j, and fij is its unit vec- 
tor. For the first order approximation, the summation is 
over all nearest-neighbor sites. Although it might over- 
estimate Ep, it can give us a clearer qualitative physical 
picture. The E^j and Ehy terms lead to the "ice rule" 
and decide the main frame of hydrogen-bond network. 
The Ep term due to dipole moments is much smaller 
than the other two terms but has an effect on the ori- 
entation of water molecules. At the same basis vectors 
(assuming |oi| = |c?2 1 = To), different isomers with differ- 
ent orientations of water molecules have the same E^j 
and Ehy but different Ep, which makes us compare the 
Ep term only. A pair of perpendicular water molecules 
are pg/47re rg lower in energy than a pair of parallel ones, 
so the orientations of two nearest water molecules prefer 
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the perpendicular type. 

Now, we discuss why only the ice nanotubes with odd 
number of side faces possess ferroelectricity. We first clas- 
sify a local water network to four types if we just consider 
the nearest pairs of water molecules. The upper panel 
of Fig. 4 shows an example of the side faces of a (n, 0) 
ice nanotube. The type (a) has the lowest energy with- 
out polarization and is preferred when an ice nanotube is 
formed. The type (b) is polarized along the axis although 
in a small higher energy. The potential energy differences 
between the four types are very small. As an example of 
(4, 0) ice nanotube, the type (b) is only 4.8 x 10~ 2 kcal 
moK 1 (2.1 x 1CP 3 cV) per water molecule higher in en- 
ergy than type (a). The lower panel of Fig. 4 illustrates 
why a (5,0) ice nanotube is ferroelectric while a (4,0) 
ice nanotube is antiferroelectric. The reason is that a 
(5, 0) ice nanotube has one ferroelectric side face of type 
(b) while a (4, 0) ice nanotube has only antiferroelectric 
side faces of type(a). This consideration can be easily 
extended to other ice nanotubes. The origin of this type 
of ferroelectricity is the Q1D confinement and the orien- 
tational order of hydrogen bonds. We get the conclusion 
that, an ice nanotube with odd number of side faces has 
at lease one redundant polarized water molecule at every 
layer and the "ice rule" accumulates them along the axis 
and generates a whole spontaneous electric polarization. 

It is worthy to mention that, during our MD simu- 
lations, many fast quenched ice nanotubes exhibit fer- 



roelectricity whatever odd or even number of side faces, 
which lie in small higher energies than those shown above. 
The side faces of these isomers are mainly type (b) and 
very stable at ice phase. So, if we apply an external 
electric field and cool the water fast, we might get ice 
nanotubes with larger polarization than those shown in 
this paper. It is also found that some (5, 0) and (7, 0) ice 
nanotubes with defects are slightly helix and ferroelectric 
too. In experiment, the orientations of water molecules 
can be detected by neutron scattering or infrared spectra, 
and the ferroelectricity can be reflected by second-order 
nonlinear optical coefficients. 

In summary, ice nanotubes inside SWCNTs with odd 
number of side faces can exhibit ferroelectricity. These 
nanoscale Q1D ferroelectric fiber structural materials 
are different from traditional ferroelectric bulk materials 
both in ferroelectric-paraelectric phase transition behav- 
ior and the original mechanism. They may be applied 
in many fields, such as high sensitive senors, nanoscale 
electric machines, nonvolatile memory devices, and so on. 
In addition, the similar type of ferroelectric mechanism 
probably exists in other nanoscale porous materials. 
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